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ABSTRACT 

Aims. 1156+295 is a flat-spectrum quasar that is loud at both radio and y-ray wavelengths. Previous observations of the source re- 
vealed a radio morphology on pc to kpc scales consistent with a helical jet model. We carried out VLB A observations at centimeter 
and millimeter wavelengths to study the structure of the innermost jet and understand the relation between the helical structure and 
astrophysical processes in the central engine. 

Methods. The source 1 156+295 was observed with the VLBA at 86, 43, and 15 GHz at four epochs from 10 May 2003 to 13 March 
2005. The observations were carried out in a full polarization mode. The highest resolution of the observations is 0.08 mas (~0.5 pc) 
at 86 GHz. 

Results. A core-jet structure with six jet components is identified in 1156+295. Three jet components are detected for the first time. 
The apparent transverse velocities of the six jet components derived from proper motion measurements range between 3.6 c and 
1 1.6 c, suggesting that highly relativistic jet plasma moves in the direction close to the line of sight. The overall jet has an oscillating 
morphology with multiple curvatures on pc scales, which might be indicative of a helical pattern. Models of a helical jet are discussed 
in the context of both Kelvin-Helmholtz (K-H) instability and jet precession. The K-H instability model is in closer agreement with 
the observed data. 

Conclusions. The overall radio structure on scales from sub-pc to kpc appears to be fitted with a hydrodynamic model in the funda- 
mental helical mode of Kelvin-Helmholtz (K-H) instability. This helical mode with an initial characteristic wavelength of A =0.2 pc 
is excited at the base of the jet on the scale of 0.005 pc (or 10 3 R S , the typical size of the broad line region of a supermassive black 
hole of 4.3 x 10 8 M o ). A precessing jet model can also fit the observed jet structure on scales of between 10 pc and 300 pc. However, 
with the precessing jet model, additional astrophysical processes may be required to explain the bendings of the inner jet structure (1 
to 10 pc) and re-collimation of the large-scale jet outflow (>300 pc). 

Key words, galaxies: jets - galaxies: kinematics - galaxies: quasars: individual: 1156+295 - radio continuum: galaxies 



1. Introduction 

Radio-loud AGNs have been observed for more than three 
decades, and tremendous progress has been made in understand- 
ing the underlying physical processes with detailed modeling 
(e.g. Urry & Padovani |19 95). However, many fundamental and 
important problems remain uncertain. For example, it is unclear 
how the jets are formed, accelerated, and collimated, and what 
causes the broad range of observed morphologies of jets. The de- 
tailed observations of structures, polarizations, and proper mo- 
tions of jets on parsec and sub-parsec scales will provide critical 
insight into how the observed morphologies of radio jets are re- 
lated to the dynamic processes of radio outflows and accretion 
disks around single, perhaps binary, supermassive black holes. 
The millimeter VLBI appears to be a unique tool to image the ex- 
tremely compact region of radio sources associated with AGNs 
at the resolution of sub-milliarcsecond, providing a wealth of in- 
formation about the detailed jet physics. 

The blazar 1156+295 (z = 0.729) is an archetypal radio- 
loud AGN with a prominent jet. It is variable across the entire 
spectrum. Strong y -ray flares have been det ected by EGRET 
and Fermi at energies higher than 100 MeV ( [Thompson et al. 



19951 [M ukherje e et al. |1997[ |Hartman et al. |1999| |Abdo et al. 
2009]). In the optical domain, the source has been classified as 
an optically violent var iable (OVV) source and hig hly polarized 
quasar (HPQ) ( [Wills et al. |1983][T^2t|Glassgold et al. |1983) . In 
the radio domain, 1156+295 displays a significant variability on 
timescales from years to days (e.g. see the monitoring program 
of the University of Michigan Radio Astronomy Observatorxn 
and the calibration program of the NRAO Very Large Arra>Q 
ILovell et al. |2003l|Savolainen & Kovalev |2008) . 

On a wide range of linear scales from pc to kpc, the source 
has a "core-jet" structure. The jet is curved and consists of sev- 
eral bright knots. On pc scales, the jet initially points to the north 
and then bends to the northeast at a distance of 3-4 mas from the 
core. At about 25 mas from the core, the jet then bends north- 
west and aligns itself with the initial direction of the kpc- scale 
jet. On the kpc scales, the jet has a wiggling morphology around 
a mean position angle of -18° an d terminates at a hotspot 2 arc- 
sec north of the core ( Ho ng et al. |20 04 ). On the basis of the high 
apparent super-luminal velocities (up to 15 h~ l c, where h is the 



Hubble constant in units of 100 km s Mpc , McHardy et al. 
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1990b, VLBI jet components, and ultra- high brightness tempera- 
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tures(> 10 12 Z)ofthe core component ( Piner & Kingham 1997 
Jorstad et al. |2001[|Hong et al. |2004| ), the angle between the jet 



axis and the line of sight is quite small. The complex bending 
jet can be attributed to the projected shape of a 3-D helical jet 
(Ho ng et al. |2004|), which appears to be common in radio-loud 
AGNs (see |Conway & Murphy |1993| ). 

To identify the relation between the wiggling jet and the ac- 
tivity of the central engine on sub-pc scales, we conducted 4- 
epoch VLBI observations of 1156+295 at three frequencies of 
86, 43, and 15 GHz (3 mm, 7 mm, and 2 cm, respectively). 
On the basis of observed morphologies, structural variations, 
and polarization properties, kinematic models based on different 
physical mechanisms are presented. Cosmological parameters of 
H = 71 km s _1 Mpc -1 , £l M = 0.27, and Q A = 0.73 are used 
throughout this paper. At the distance of 1 156+295, 1 mas corre- 
sponds to 7.3 pc, and a proper motion of 1 mas yr" 1 corresponds 
to 23.7 c ( [Wight |2006| ). 

2. Observations and data reduction 

The multi-frequency observations were carried out with the Very 
Long Baseline Array (VLBA) on 10 May 2003 (2003.36), 24 
July 2003 (2003.56), 1 April 2004 (2004.25), and 13 March 2005 
(2005.19) at 86, 43, and 15 GHz. Detailed information about the 
observations is given in Table [T] 

We used a snapshot observing mode (one scan of 4-5 min- 
utes duration) for our main target source and calibrators. To opti- 
mize the u-v coverage, the three observing frequencies of 86, 43, 
and 15 GHz were cycled in each of the observing runs. A point- 
ing scan was inserted after each of the frequency switching to 
86 GHz. A total observing time around 6 hours on 1156+295 
was evenly allocated to three observing frequencies at each 
epoch. 

All observations were made in the dual circular-polarization 
mode. At 15 and 43 GHz, 3C273, J13 10+3220 and OJ 287 were 
included in the observing runs as calibrators for corrections for 
R-L delay-offset, polarization leakage, and the absolute position 
angle of the electric vector (EVPA), respectively. At 86 GHz, it 
is difficult to find a compact calibrator with a stable polarization 
angle to use for EVPA calibration, hence we only study the to- 
tal intensity structure of 1156+295 from the 86 GHz data. The 
raw data were recorded in the 2-bit VLBA format with a total 
bandwidth of 64 MHz and correlated in Socorro, New Mexico, 
USA. 

The data were reduced following the standard procedure 
of VLBA data reduction using the NRAO AIP S (Astronomica l 
Image Processing System) software package (Dia mond |1995| ). 
The visibility amplitudes were calibrated using the gain curves 
and system temperature monitored during the observations. The 
atmospheric opacity correction was applied by using the WX ta- 
ble, which is an extension file of data in a table format contain- 
ing ground weather information. Fringe fitting was made using a 
short scan on 3 C 273 to calibrate the delays and phases. The so- 
lutions were then interpolated to the data on all the programme 
sources. Global fringe fitting was carried out to help correct for 
the residual delays and fringe rates. At 15 and 43 GHz, 3C 273 
and OJ 287 were used to calibrate the complex bandpasses. 

A series of polarization-specific calibrations were made for 
the 15 and 43 GHz data. After the amplitude calibration, vari- 
ations in the parallactic angles were determined and removed 
from the data. The R-L delay difference was determined by run- 
ning the VLBA procedure VLBACPOL on the highly polarized 
source 3C 273. The instrumental polarization parameters of the 



antenna feeds were calibrated with the AIPS task LPCAL us- 
ing scans on J13 10+3220. The corrections were written in an 
updated antenna (AN) table and then applied to the data. The 
calibration of the absolute polarization angle was made from 
observations of OJ 287 or J13 10+3220. At 15 GHz, we used 
OJ 287, which is compact and stable enough to be an EVPA cal- 
ibrator. To obtain the EVPA correction, we compared its appar- 
ent EVPA with the near- simultaneous measurements included in 
the MOJAVE/2cm Survey databas^] At 43 GHz, both OJ 287 
and Jl 3 10+3220 were used for similar calibrations. Corrections 
were derived by comparing their apparent EVPA with the 
VL A/VLB A measurement^] at 43 GHz on close epochs. After 
calibration, the data were divided into single-source file s and im- 
ported into the Caltech VLBI Program DIFMAP ( [Shepherd etaL 



1994 ) for self-calibration and imaging. The fully calibrated data 
at 15 and 43 GHz were then re-imported into AIPS to generate 
polarization images. I, Q, and U images were produced sepa- 
rately with the AIPS task IMAGR using natural weighting. The 
polarization angle and polarization intensity images were then 
obtained by combining the Q and U images using the AIPS task 
COMB. 

At 86 GHz, an antenna-based edition of the data was first 
produced in DIFMAP to remove visibilities related to extremely 
large amplitude fluctuations likely due to bad weather or imper- 
fect atmospheric opa city corrections. Th e uv data were then av- 
eraged using Miriad ( [Sault et al~||1995| ) over a time span of 1 
minute. We calculated the scalar average in amplitude because 
the 86-GHz visibility phase is significantly affected by fast atmo- 
spheric path length changes. After that, the non-imaging analysis 
was done in the wv-domain (see the Section [3X5]). Finally, the 
86-GHz data were re-imported into DIFMAP for self-calibration 
and imaging. The image at 86-GHz was made with uniform 
weighting to achieve a higher angular resolution. A few iter- 
ations of phase-only self-calibration were made until the rms 
noise in the image did not decrease. Amplitude self-calibration 
was not performed because we found that it would cause an un- 
certainty as large as 20-30 percent in the flux density at 86 GHz. 



3. Results and analysis 

3.1. Structure of the source 1 156+295 
3.1.1. Images at 15 and 43 GHz 

The contour images in Figures [T] and [2] represent the distribution 
of total radio intensity in 1156+295 at four epochs at 15 and 
43 GHz, respectively. The parameters of the images are listed in 
Tabled 

In the images at 15 GHz, the peak values of the intensity 
are 1.61, 1.86, 0.46, and 1.31 Jy beam -1 for the four epochs, 
respectively. At 43GHz, the values are 1.54, 2.22, 0.44, and 
1.59 Jy beam -1 , respectively. The peak value at the highest 
epoch (2003.56) is four or five times as high as at the lowest 
epoch (2004.25). This large variation in peak intensity is not 
due to calibration error, but reflects an intrinsic variation in the 
source. The similar trend in the flux density variation is seen in 
the flux density monitoring program at centimeter and millimeter 
wavelengths {e.g. VLA: 5, 8.4, 22, and 43 GHz, Motsahovi: 22 
and 37 GHz, OVRO MMA: 1 10 GHz), confirming that the source 
started flaring around epoch 2003.3, reached the peak around 
epoch 2003.6 and at higher frequencies earlier, and that the total 



3 http://www.physics.purdue.edU/MOJAVE/sourcepages/l 156+295. shtml 

4 http ://w w w. vla.nrao . edu/astro/calib/polar 
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Obs. date 


Epoch 


Frequency 


VLBA 


Band Width 


On- source Time 




(yr) 


(GHz) 


participated telescopes 


(MHz) 


(Min) 


2003-05-10 


2003.36 


15.0 


All 10 


64 


122 






43.0 


9 of 10 a 


64 


122 






86.0 


7 of 10 b >° 4 


64 


121 


2003-07-24 


2003.56 


15.0 


9 of 10 e 


64 


122 






/ion 
43. U 


9 OI 10 


54 


1 n 

Izz 






86.0 


8 of 10 ^ 


64 


121 


2004-04-01 


2004.25 


15.0 


All 10 


64 


118 






43.0 


All 10 


64 


118 






86.0 


8 of 10 b > c 


64 


117 


2005-03-13 


2005.19 


15.0 


9 of 10 1 


64 


118 






43.0 


8 of 10 ^ 


64 


118 






86.0 


8 od 10 °J 


64 


117 



Notes: a: FD data were corrupted at 43.0 GHz; b: BR had no 3 mm receiver; c: SC had no 3 mm receiver; d: No fringes detected at 86.0 GHz on 
baseline of HN; e: SC was out for maintenance; /: HN was removed during the observing run due to bad weather; g: SC data were flagged out 
because of their abnormal system temperature. 



flux density decayed rapidly to its minimum value around epoch 
2004.4. The next flare appeared to reach its peak around 2005.4. 

In 15 GHz images, the source has a compact and bright 
region (referred to as the "core" hereafter) that dominates the 
brightness distribution and has a one-sided jet to the north. The 
jet initially extends to the northeast, bends to the northwest and 
becomes a diffuse structure at a distance of 3-4 mas from the 
core. At about 12 mas from the core, the brightness of the jet 
drops below the detection limit of 0.6 mJy beam _1 (3<x). The 
overall structure of 1156+295 observed at 15 GHz is in good 
agreement with the previous observations at lower frequencies 



( |Hong et al. |200^[Pmer & Kingham 1 997l|Jorstad etal. |2001| ). 

At 43 GHz, with an angular resolution of 0.2 mas (corre- 
sponding to 1.5 pc, and close to the size of the broad line region 
of a normal AGN), our observations reveal a core-jet structure 
on a scale smaller than 1 mas (7.3 pc). The jet first extends to 
the north in the innermost 0.5 mas, then bends to the north-east 
along the direction at PA ~18° with respect to the core. The mor- 
phology is consistent with the structure observed at 15 GHz. 



188.5 mJy. In comparison to C5, both the size and flux density 
of C6 varied relatively slowly as it moved out from the core. 

At 43 GHz, in the first three epochs (2003.36, 2003.56, and 
2004.25), both C5 and C6 were detected. In addition, C7, the 
closest component to the core was detected (see the right panel of 
Figure [3}. All the components C5, C6, and C7 were detected for 
the first time at the epoch of 2003.36 at this frequency. We found 
C6 varied significantly and the variation in the flux density at 
43 GHz appeared to be correlated with that observed at 15 GHz. 
We also note that from 2003.36 to 2003.56, the apparent flux 
density of C7 increased by a factor of 2.3 and the apparent size 
increased by a factor of 2.7. The possible reason for the rapid 
increase in the radiation intensity of C7 is the absence of the SC 
(San Croix) antenna which decreases the resolution and causes 
an increase in the apparent size, hence the apparent flux density 
of C7. In the fourth epoch (2005.19), only two components were 
detected because of the relatively poorer sensitivity and angular 
resolution in the absence of two VLBA antennas (HN and SC); 
the data were fitted with two Gaussian components. 



3.1 .2. Model fitting of data at 15 and 43GHz 



3.1 .3. The visibility analysis and images of 86GHz data 



To quantitatively analyze the properties of 1156+295, we fitted 
the source emission structure with Gaussian components using 
the program MODELFIT in DIFMAP. The model-fitting param- 
eters are listed in the left panels of Tables [3] and [4] for 15 and 
43 GHz, respectively. 

At 15 GHz, in addition to four sets of 15 GHz data described 
in Section [2j three additional data sets for epochs 2001.17, 
2002.90, and 2004.11 from the "MOJAVE 2cm Survey Data 
Archive" are used in the model fitting. The source structure is 
most closely fitted with six Gaussian components (see the left 
panel of Figure [3]) at this frequency. The brightest region (com- 
ponent C) could be fitted with the elliptical Gaussian component 



that was considered to be the core in the previous works ( Hong 
et al. |2004| ). The jet is fitted with five circular Gaussian compo 



nents, which are named 'C2' to 'C6' following the nomenclature 
in Hong et al. (2004 ). The component C5 was first detected at 
the epoch 2001.17 with a flux density of 1 15.3 mJy at a distance 
of ~0.7 mas in PA~13°. Its flux density dropped quickly to about 
10% of its initial value at the epoch 2005.19 and the size of the 
component expanded slightly as C5 moved out from the core. 
C6 was first detected at the epoch 2003.36 with a flux density of 



We inspected the visibility amplitude of 86 GHz data as a func- 
tion of the projected baseline to acquire a general idea of the 
source structure. The visibility amplitude profiles are not in con- 
sistent with the expectations for a simple point source, but there 
is some evidence of intrinsic source structures. 

We fitted the amplitude profiles on three different baseline 
length ranges, the shortest (60-70x1 6 /i), intermediate (100- 
200xl0 6 ^), and longest (1100-1500xl0 6 ^), with a point source 
model, respectively. The correlated flux density on the short- 
est baselines (S s ) is 1.5, 1.6, 0.8, and 0.7 Jy at epoch 2003.36, 
2003.56, 2004.25, and 2005.19, respectively (see Table[5]), sug- 
gesting a flux density variation on a timescale of years to months 
(the shortest time separation being 2.4 months). The largest vari- 
ation accounts for a factor of 2.3. The trend in the flux density 
variation observed at 86 GHz agrees with that seen for the 15 and 
43 GHz observations at these epochs, and the variations in flux 
density reported from the monitoring program at centimeter and 
millimeter wavelengths (see Section 3.1.1| ). The correlated flux 
densities on the intermediate baseline (S m ) and longest baseline 
(Si) are lower than S s , providing evidence that the source struc- 
ture at 86 GHz is resolved by the VLBA at these epochs. 



4 W. Zhao: Sub-pc structure of blazar 1 156+295 



Table 2. The parameters of contour images at 15 and 43 GHz in Figures 



Frequency 


Epoch 


Synthetized Beam 


-^peak 


lcr 


Contours 








Maj 


Min 


P.A. 










(GHz) 


(yr) 


(mas) 


(mas) 


(°) 


(Jy beam -1 ) 


(mJy beam -1 ) 


(mJy beam -1 ) 




15.0 


2003.36 


0.95 


0.57 


-9.5 


1.61 


0.20 


0.60 x (-1,1,2,4, 


.., 1024) 




2003.56 


1.34 


0.74 


-21.7 


1.86 


0.30 


0.90 x (-1,1,2,4, 


.., 1024) 




2004.25 


0.94 


0.58 


-11.0 


0.46 


0.20 


0.60 x (-1,1,2,4, 


.., 1024) 




2005.19 


1.17 


0.67 


13.0 


1.31 


0.20 


0.60 x (-1,1,2,4, 


1024) 


43.0 


2003.36 


0.35 


0.19 


-11.8 


1.54 


0.40 


1.20 x (-1,1,2,4, 


1024) 




2003.56 


0.47 


0.27 


-23.2 


2.22 


0.60 


1.80 x (-1,1,2,4, 


.., 1024) 




2004.25 


0.33 


0.20 


6.2 


0.44 


0.30 


0.90 x (-1,1,2,4, 


...,256) 




2005.19 


0.53 


0.34 


-10.9 


1.59 


0.80 


2.4 X (-1,1, 2, 4,.. 


,512) 



Table 3. 15 GHz model-fitting results and physical parameters of fitted models 



Epoch 


Comp 


S c 


r 


PA 


a 


b 1 a 


O 


T h 


<5eq 




y 


6 


(yr) 




(mJy) 


(mas) 


(°) 


(mas) 




(°) 


K 




K 




O 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


2001.17 


C 


996.2+149.4 








0.21 


0.16 


43.4 


2.0 x 10 12 


8.8 


2.2 x 10 11 








C5 


115.3 + 17.3 


0.69+0.03 


13 +3 


0.35+0.03 


1 





1.3 x 10 10 






5.2 


4.5 




C4 


6.2 + 1.6 


1.31+0.62 


-14+27 


0.22+0.62 


1 





1.8 x 10 9 






5.5 


5.3 




C3 


30.7 + 32.0 


4.32+ 1.91 


-8 + 24 


3.67+3.81 


1 





3.1 x 10 7 






- 


- 




C2 


15.9 + 14.7 


6.69+0.89 


18 + 8 


1.96+1.79 


1 





5.7 x 10 7 






- 


- 


2002.89 


C 


2065.5+309.8 








0.37 


0.31 


4.7 


0.7 x 10 12 


2.3 


2.9 x 10 11 








C5 


48.5 + 7.6 


0.69+0.11 


23 +9 


0.23+0.11 


1 





1.2 x 10 9 






4.2 


20.5 




C4 


9.5 + 2.6 


1.86+0.58 


-16+18 


0.52+0.58 


1 





4.8 x 10 8 






5.4 


19.1 




C3 


35.9 +25.5 


4.51+0.87 


9+11 


2.48+1.74 


1 





7.9 x 10 7 






- 


- 




C2 


31.6 +25.8 


8.67+ 1.44 


24 + 10 


3.57+2.88 


1 





3.4 x 10 7 








- 


2003.36 


C 


1605.5+240.8 








0.23 


0.31 


2.8 


1.3 x 10 12 


5.4 


2.5 xlO 11 








C6 


188.5 + 28.3 


0.56+0.02 


13 +2 


0.33+0.02 


1 





2.3 x 10 10 






4.7 


10.1 




C5 


41.3 + 6.2 


0.94+0.07 


30 +4 


0.33+0.07 


1 





5.1 x 10 9 






4.0 


9.5 




C4 


9.4 + 1.5 


1.98+0.31 


-12+9 


0.40+0.31 


1 





8.1 x 10 8 






4.4 


10.0 




C3 


28.1 +5.5 


4.56+0.22 


1 +3 


2.78+0.45 


1 





4.9 x 10 7 






- 






C2 


28.1 + 14.5 


8.19+1.01 


26 + 7 


3.93+2.02 


1 





2.5 x 10 7 






- 




2003.56 


C 


1907.5+286.1 








0.28 


0.31 


15.7 


1.1 x 10 12 


4.3 


2.5 x 10 11 








C6 


135.8 + 20.4 


0.59+0.04 


17 +4 


0.39+0.04 


1 





1.2 x 10 10 






4.6 


12.3 




C5 


38.7 + 5.9 


1.21+0.14 


26 +7 


0.40+0.14 


1 





3.2 x 10 9 






3.8 


12.5 




C4 


29.0 + 4.5 


2.25+0.18 


-10+5 


0.94+0.18 


1 





4.5 x 10 8 






4.4 


12.8 




C3 


28.2 + 29.1 


4.64+1.74 


2+21 


3.38+3.49 


1 





3.3 x 10 7 












C2 


24.3 + 34.4 


8.80+3.10 


27 +19 


4.38+6.19 


1 





1.7 x 10 8 










2004.11 


C 


625.7 + 93.9 








0.31 


0.23 


9.1 


0.4 x 10 12 


1.3 


2.9 x 10 11 








C6 


122.6 + 18.4 


0.57+0.02 


19 +2 


0.28+0.02 


1 





2.1 x 10 10 






8.8 


21.4 




C5 


27.2 + 4.2 


1.11+0.07 


21 +4 


0.41+0.07 


1 





2.2 x 10 9 






6.0 


24.7 




C4 


11.5 + 3.1 


1.96+0.17 


-12+5 


1.33+0.17 


1 





8.8 x 10 7 






8.1 


22.1 




C3 


21.8 + 17.8 


4.73+1.10 


2+13 


2.69+2.19 


1 





4.1 x 10 7 












C2 


31.9 + 19.1 


8.97+1.54 


29+10 


5.17+3.09 


1 





1.6 x 10 6 










2004.25 


C 


501.0 + 75.2 








0.43 


0.34 


9.5 


0.1 x 10 12 


0.3 










C6 


120.3 + 18.1 


0.62+0.01 


18 +1 


0.35+0.01 


1 





1.3 xlO 10 












C5 


22.6 + 3.4 


1.25+0.07 


22 +3 


0.50+0.07 


1 





1.2 x 10 9 












C4 


6.3 + 1.1 


2.33+0.23 


-18+6 


0.60+0.23 


1 





2.4 x 10 8 












C3 


26.0 + 19.4 


4.54+1.18 


0+15 


3.17+2.36 


1 





3.5 x 10 7 












C2 


25.3 + 13.3 


8.94+1.44 


32 + 9 


5.5 +2.88 


1 





1.1 x 10 7 










2005.19 


C 


1307.9+196.2 








0.08 


1 





2.5 x 10 12 


11.4 


2.2 x 10 11 








C6 


64.5 + 9.7 


0.89+0.02 


15 +1 


0.33+0.02 


1 





7.9 x 10 9 






6.6 


3.4 




C5 


11.9 + 1.8 


1.27+0.09 


20 +4 


0.40+0.09 


1 





1.0 x 10 9 






6.3 


2.9 




C4 


5.0 + 1.1 


2.43+0.43 


-11 + 10 


1.16+0.43 


1 





5.1 x 10 7 






6.6 


3.3 




C3 


26.1 + 8.5 


4.09+0.33 


4+5 


2.07+0.67 


1 





8.3 x 10 7 












C2 


48.8 + 11.3 


9.43+0.79 


28 +5 


6.94+1.58 


1 





1.4 x 10 7 











Note: S c is the flux density of the fitted components; r and PA are the radial separation and the position angle of the jet components with respect to 
the core component, respectively; a and b are the major and minor axis size for the Gaussian model in mas; T h is the brightness temperature in the 
parent-galaxy rest-frame, £ eq is the equipartition Doppler factor and T r is the intrinsic brightness temperature eliminating the Doppler boosting; y 
and 6 are the Lorentz factor and the viewing angle of jet components. 
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Table 4. 43 GHz model-fitting results, and physical parameters of fitted models 
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Epoch 


Comp 




r 


PA 


a 


b/a 


o 


T h 






y 


6 


(yr) 




(mJy) 


(mas) 


(°) 


(mas) 




(°) 


K 




K 




(°) 


(1) 


U) 






(?) 


w 


\l) 


\P) 


\y) 


(10) 


(11) 


(12) 


(1J) 


2003.36 


c 


1539.5+230.9 








0.05 


0.15 


-30.3 


6.8 x 10 


13.4 


5.1 X 10 11 








C7 


75.5 + 11.4 


0.22 


15 


0.06 


1.00 





3.5 x 10 






7.3 


2.2 




C6 


105.3 + 15.8 


0.46 


4 


0.30 


1.00 





1.9 x 10 9 














41.4 ± O.J 


n en 


Ad 


n 9c 


1 nn 


u 


o o w i r»8 
o.o X 1U 










2003.56 


C 


2177.5+326.6 








0.05 


0.42 


-34.2 


3.4 x 10 


16.5 


2.1 x 10 








CI 

/ 


177 zl -i- 9^ Q 
1 /Z.4 ± zj.y 


n 9Q 


Zl 


n 1 ^ 

U. 10 


1 nn 
1 .uu 


n 
u 


1 1 v 1 ni° 

1.1 X 1U 






o. / 


i .j 




C6 


62.9 + 9.6 


0.63 


14 


0.44 


1.00 





5.4 x 10 8 












C5 


9.9 + 2.1 


1.11+0.03 


26 +2 


0.41+0.03 


1.00 





2.9 x 10 7 










2004.25 


C 


482.6 +72.4 








0.10 


0.40 


-1.2 


0.2 x 10 12 


0.6 










C7 


38.4 + 5.8 


0.38 


-3 


0.12 


1.00 





4.4 x 10 9 












C6 


61.2 + 9.2 


0.62 


18 


0.23 


1.00 





1.9 x 10 9 












C5 


27.5 + 4.3 


1.17+0.01 


16+1 


0.40+0.01 


1.00 





2.9 x 10 8 










2005.19 


C 


1648.3+247.3 








0.10 


0.60 


-89.0 


0.5 x 10 12 


6.8 


0.7 x 10 11 








C6 


24.4 + 4.1 


0.78+0.02 


12+2 


0.27+0.02 


1.00 





5.5 x 10 8 











Note: S c is the flux density of the fitted components; r and PA are the radial separation and the position angle of the jet components with respect to 
the core component respectively; a and b are the major and minor axis size for the Gaussian model in mas; T h is the brightness temperature in the 
parent-galaxy rest-frame, £ eq is the equipartition Doppler factor and T r is the intrinsic brightness temperature eliminating the Doppler boosting; y 
and 6 are the Lorentz factor and the viewing angle of jet components. 



We then use the brightness distribution models introduced by 
Pearson (1995|) to estimate the flux density (S /) and size (6^) of 



the source. For the first three epochs, an extended source model 
can be used (see Figure [4]). The calculation gives the flux den- 
sity 1.7, 1.8, and 0.8 Jy, respectively, and the largest size of the 
source 0.06 mas (2003.36), 0.05 mas (2003.56), and 0.05 mas 
(2004.25), respectively (see Table [5]). For the first three epochs, 
the correlated flux density on the shortest baseline (S s ) is com- 
parable to the fitted Gaussian peak of the amplitude profile, im- 
plying that S s is a reasonable representative of the zero-spacing 
flux density. For the last epoch (2005.19), a double-point source 
model was used (see Figure]?]). The flux density is 0.9 Jy, the 
flux density ratio of two point sources is 3:1, and the separation 
of the two point components is 0.16 mas (see Table [5]). 

On 86 GHz CLEANed images (see Figure [5]), most of the 
jet emission is resolved out at an angular resolution of 0.2 mas 
xO.08 mas. At the epoch 2003.36, a weak jet is seen ~0.3 mas 
north of the core. The 86 GHz jet coincides with the inner- 
most jet component C7 detected at 43 GHz at the same epoch. 
Earlier work, based on Global Millimeter VLBI Array observa- 
tions on 26 October 2001 at 86 GHz with an angular resolution 
of 0.188x0.038 mas, detected a jet component at 0.19 mas north 
of the core, marginally coinciding with C7 in our observations. 
At all other three epochs of 2003.56, 2004.25, and 2005.19, the 
jet emission on scales of between 0.1 and 1 mas detected at 
both 43 and 15 GHz is below the detection threshold at 86 GHz 
(i.e., 3<x 30m Jy beam' 1 ), showing only a point-like source in 
1156+295. 



3.2. Flux variability of the core component 

The radio emission of the source is dominated by the compact 
core at all three observing frequencies of our experiment. The 
variability of the core component C at each frequency can be 
assessed using the index of variability V 



Table 5. 86 GHz parameters from model fitting of amplitude 
profiles 



Epoch 


s s 




Si 


Sf 


o d 


T h 


(yr) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


K 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


2003.36 


1.5 


1.1 


0.6 


1.7 


0.06 


1.9 x 10 11 


2003.56 


1.6 


1.2 


1.3 


1.8 


0.05 


2.9 x 10 11 


2004.25 


0.8 


0.4 


0.3 


0.8 


0.05 


1.2 x 10 11 


2005.19 


0.7 


0.5 


0.7 


0.9 a 


0.16 




Errors 


20% 





Note: a it is fitted with a double point- source model, the separation of 
the two points is 0.16 mas, the flux density ratio is 3:1. 
S 5 , S m , and Si are the correlated flux density on the short- 
est (60-70xl0 6 ^), intermediate ( 1 00-200 xlO 6 A), and longest (1100- 
1500xl0 6 /i) baselines, respectively; S f and 6 d are the flux density and 
the size of the fitted Gaussian model; T h is the brightness temperature 
in the parent-galaxy rest-frame. 



where S max and S m [ n are the maximum and minimum fitted flux 
densities of the core component. A value of V ~ represents no 
variability and V ~ 1 extreme variability. During the observing 
period, V was 0.58+0.07, 0.64+0.06, and 0.39+0.13 for 15, 43, 
and 86 GHz, respectively, suggesting that the component C is 
highly variable (the uncertainties in V are evaluated based on 
the uncertainties in the S max and S m i n ). 

The core of an AGN is often defined to be the most com- 
pact region of the jet or the jet base, which always appears 
to be optically thick at centimeter wavelengths because of 
the superposition of sub-sections of the inner jet with various 
synchrotron self- absorption turnover frequencies. At millime- 
ter wavelengths, many quasars and AGNs display a spectral 
turnover at 10 11±0 5 GHz (32-320 GHz)(|Marscher |1995|). In this 



V = 



* max $ YYiin 

* max ~l~ S min 



(1) 



frequency range, the inner jet becomes optically thin, allowing 
us to observe the intrinsic source structure with VLBI at mm- 
submm wavelengths. Figure [6] shows the flux density of the core 
at three frequencies at each of four epochs. It shows that at the 
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epochs 2003.56 and 2005.19, the flux density at 43 GHz dom- 
inates over those at 15 and 86 GHz, implying that the turnover 
frequency is around 43 GHz. 



3.3. Superluminal motion of jet components 

The angular distances to the core of all detected jet components 
were measured at all observing epochs. We carried out linear re- 
gression fitting of the separation of the jet components from the 
core as a function of the observing epoch. From the fitted slope 
of the linear regression, we determined the apparent velocities of 
each jet component. 

The uncertainty in the separation of jet components to the 
core component C is estimated to be AR = (for the compact 

components from C4 to C7), or AR = (for extended 

components C2 and C3), where W is the FWHM of the synthe- 
sized beam, a is the major-axis size of the componen t, and SNR 
is the signal-to-noise ratio for a given component ( [FomalonF 
|2004l|Lee et al. |2008| ). 



To better determine the proper motions of the components 
C2, C3, and C4, in addition to the data discussed in the present 
paper, we also included pub lished data at low er observing fre- 
quencies at the early epochs ( Hon getal. 2004| ). For components 
C5 and C6, both 15 and 43 GHz data from this work were used 
in the analysis. Since C7 can be separated from the core emis- 
sion at only 43 GHz, the determination of the proper motion for 
C7 involves only the 43 GHz data from this paper. 

We use the inverse square of the uncertainty, (l/(Ar) 2 ), as a 
weight parameter in the least squares fitting for C2, C3, C4, and 
C7. For C5 and C6, we used both 15 and 43 GHz data points 
in the least squares fitting. The inverse square of the difference 
between the fitted positions at 15 and 43 GHz was used as a 
weight parameter. 

Figures [7] and [8] illustrate the linear fit of the angular separa- 
tion of the jet components for the determination of their proper 
motions. Significant proper motions of 0.49+0.03, 0.18+0.03, 
0.15+0.02, 0.19+0.02, and 0.16+0.04 mas yr" 1 were determined 
for the jet components C2, C4, C5, C6, and C7, respectively. 
The apparent transverse velocities of 11.6+0.7 c, 4.3+0.7 c, 
3.6+0.5 c, 4.5+0.5 c, and 3.8+0.9 c derived from the proper mo- 
tions for C2, C4, C5, C6, and C7, respectively, confirming early 
results on a pparent superluminal motion in the blazar 1 156+295 
(Hong et al. |[2004| . We note that C3 shows a peculiar pattern 
in its proper motion. The fitting to the data before 2001 gives a 
proper motion of 0.41+0.08 mas yr" 1 (or 9. 7+ 1.9c). No signifi- 
cant proper motions was detected after 2001.17 when the com- 
ponent reached the place where the jet changes its direction by 
~ 40°. Both the proper motion velocities and the component 
sizes of the jet components show an increase beyond a distance 
>2 mas (or 14 pc) away from the core. A possible explanation of 
these changes is a decrease in the external pressure gradient. 

We also note that the inferred values of apparent speed for 
C4, C5, C6, and C7 are similar, and the apparent speeds of C2 
and C3 are much higher than those of C4 to C7. A possible ex- 
planation is that the extended components C2 and C3 are rather 
diffuse, and the surface brightness of these components is very 
low. The uncertainty in the radial distance from the core is large 
and might affect the fitting value of the apparent speed. 



3.4. Physical properties of compact components in VLBI 
images 

The brightness temperature in the parent-galaxy rest-frame 7\> 
can be derived from the measurements of the flux density (Sob), 
the observing frequency (v b), and the angular size (6$) after cor- 
rection for the cosmological expansion effect (1 + z) ( Giiijosa & 
|Daly 1 1996) ) 



1.77 x 10 12 (1 +z) 



>ob 



"ob' 



(2) 



where Sob * s m Jy» y ob in GHz, and 6d is in mas. In our paper, 
Qd = ^fab for a Gaussian component with major (a) and minor 
(b) axis sizes at 15 and 43 GHz. 

The derived results of 7b are summarized in the Col. 9 of 
Tables [3] and |4| and Col. 7 of Table|5] For some epochs at 15 and 
43 GHz, we note that the observed brightness temperature 7b 
of the core component is close to or exceeds the upper limit of 
7 b restricted by the "inverse Compton cata strophe" ( Keller man & 



Pauliny-Toth 1969 ) or equipartition limit ( Readhead 1994 ). The 



higher brightness temperature > 10 12 K in flat-spectrum radio 
sources is often considered to be a result of the Doppler boosting 
effect and /or indications on a non- stationary process. 

We note that at the epochs 2003.36 and 2003.56, the bright- 
ness temperatures at 86 GHz are much lower than that mea- 
sured at 43 GHz. This phenomenon seems to be common in 
flat-spectrum radio-loud quasars and AGNs ( Lee et af~||2008| 
Lobanov et al. | [2000). Since 1156+295 is still rather compact 



at 86 GHz, e.g. the ratio of the CLEANed flux density to the 
correlated flux density is approximately equal to 1 in our ob- 
servations. Therefore, the relativel y lower 7\, at 86 GHz is not 
caused by the loss of compactness ( |Lee et al."|20 08) but may be 
caused by the intrinsic changes in the physical conditions of the 
innermost jet (<0.04 mas, or 0.30 pc), e.g. the emission structure 
becomes optically thin at 86 GHz and the brightness temperature 
is intrinsically lower than the optically thick case. 

To estimate the Doppler factor of 1156+295 using the re- 
sults of our observations, we assume the equipartition of the en- 
ergy between the radiating particles and the magnetic field. The 
equipartition Doppler factor of a uniform self-absorb ed source 
can be obtained from Eq. 4 in |Guijosa & Daly | ( |1996| ) if we as- 
sume that the emission arise from an optical-thin, homogeneous 
sphere. 

£ eq = [[10 3 F(a)] 34 {[l - (1 + zr l/2 ]/(2h)}- 2 S lp- 34 x (1 + 

z )15-2or( y x 1 Q3^-(2flr+35)jl/(13-2or) 

where h is the Hubble constant in units of 100 km s" 1 Mpc -1 , 
the spectral index a is assumed to be -0.75, and 0d is the an- 
gular diameter of the homogenous sphere in mas. For elliptical 
Gaussians, 0d = 1.8 ^fab, where a is the major axis size and b 
is the minor axis size. In addition, v op is the turnover frequency 
in GHz. To simplify the assessment, we used the observing fre- 
quency 15 and 43 GHz for replacing the turnover frequency. The 
data at 86 GHz are not used here because the source might be 
transparent to the synchrotron emission at this frequency. Thus, 
a higher limit to the equipartition Doppler factor is inferred from 
our estimate, where S op is the flux density in Jy at the turnover 
frequency. In this paper, we use the observed flux density at 15 
and 43 GHz as the S op . F(a) is a factor with a value approxi- 
mately equals to 3.4 at the turnover frequency. 

The results are listed in column 10 of Tables [3] and 0] The 
values of the inferred Doppler factor are in the range of 1.3 to 
16.5 except for the epoch 2004.25 when the source is in the low- 
est state. 



W. Zhao: Sub-pc structure of blazar 1156+295 
Table 6. The parameters of contour images in Figured 
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Epoch 


Synthetical Beam 


^peak 


lcr 


Contours 




Sc 




Maj 


Min 


P.A. 












(yr) 


(mas) 


(mas) 


(°) 


(Jy beam -1 ) 


(mJy beam -1 ) 


(mJy beam -1 ) 




(Jy) 


2003.36 


0.22 


0.08 


-24.3 


0.92 


4.99 


15.0 X (-1,1, 2, 4,.. 


.,32) 


1.27 


2003.56 


0.23 


0.08 


-24.7 


1.96 


10.00 


30.0 x (-1,1, 2, 4,.. 


.,64) 


2.02 


2004.25 


0.21 


0.08 


-25.1 


0.58 


4.20 


12.6 x (-1,1, 2, 4, 8 


,16) 


0.58 


2005.19 


0.24 


0.11 


-11.4 


0.73 


5.90 


17.7 x (-1,1, 2, 4,.. 


.,64) 


0.74 



Note: S c is the CLEANed flux density 



Table 7. Linear polarization properties at 15 and 43 GHz 



Frequency 


Epoch 


Component 


•^peak 


F P 


m 


EVPA 


(GHz) 


(yr) 




(mJy beam -1 ) 


(mJy) 


(%) 


(°) 


15.0 


2003.36 


C 


28.9+0.3 


30.4+0.7 


2.0+0.4 


112+7 






C5 


9.9+0.3 


3.7+0.2 


6.0 +2.2 


75+6 




2003.56 


C 


34.6+0.2 


36.7 +0.8 


2.2 +0.7 


106+7 






C5 


9.6+0.2 


4.4+0.2 


6.7+2.6 


74+4 




2004.25 


C 


12.8+0.1 


14.3 +0.3 


3.6+1.1 


106+2 






C5 


10.6+0.1 


5.5 +0.2 


6.7+2.0 


80+3 




2005.19 




18.2+0.3 


20.4 +0.2 


4.8+4.4 


44+3 


43.0 


2003.36 




36.6+0.5 


27 +0.8 


2.4 +0.3 


120+ 8 




2003.56 




29.1 + 1.0 


17.5 +1.0 


2.7 +0.3 


35+12 




2004.25 




5.2+0.4 


1.5 +0.2 


1.8+0.5 


72+5 



Note: 

Fpeak^ Peak polarized brightness. 
F P : Integrated polarized flux density, 
m: the polarization percentage. 
EVPA: Electric-Vector Position Angle. 



Since S Qq strongly depends on the observed values of S 



op, 

y op , and Gj, the difference between the values of S Qq at different 
frequencies for a given epoch are dominated by the error in the 
observed quantities. The intrinsic brightness temperature elimi- 
nating the Doppler boosting in the source rest- frame is T Y , and 
T v = T\j/6. Assuming 6 = £ eq , we can estimate the value of T r . 

The inferred physical quantities of T\» S tq ^ and T v are sum- 
marized in the right panels of Tables [3] |4j and|5] 

The kinematic properties (the Lorentz factor y, the viewing 
angle 6, and the intrinsic velocity j3 in the unit of c) for the com- 
pact components C4 to C7 (the components C2 and C3 appear 
too extended to have meaningful assessments for the kinematic 
properties) in the inner j et region are dis cussed on the basis of 
Equations (3-5) given in Ghisell ini et al. (1993]) 



25 



6 = arctan ( 



2fi, 



app 



app + ^ 2 



^app 



sin 6 + y^app cos 6 ' 



(3) 



(4) 



(5) 



where J3 app is the apparent proper motio n ve locity of jet compo- 
nents derived by linear fitting in Section 3.3 and S is the Doppler 
factor with a value derived from the equipartition calculation £ eq . 

The inferred Lorentz factor ranges between 3.8 and 8.8, the 
viewing angle ranges between 1.5° and 24.7° (see Tableland 



[4]), and the intrinsic jet speed ft approximately equals to 0.98 c. 
The jet properties inferred from the kinematic quantities confirm 
that 1156+295 is a source affected by the extreme relativistic- 
beaming effect. 



3.5. Linear polarization properties of 1156+295 

Figures [T] and [2] show the linear polarization structure of 
1 156+295 at 15 and 43 GHz, superimposed on the contour maps 
of total intensity at the corresponding frequencies. The polar- 
ization properties of the source, the polarized flux density F po i, 
peak polarization intensity Fpeak, the fractional polarization m, 
and average EVPA are given in Table [7] 

At 15 GHz, for the first three epochs (2003.36, 2003.56, and 
2004.25) there is consistent structure in polarized emission. In 
the inner part of the structure, we probably see the polarized 
emission located in the areas of the core C, component C6 and 
C7. In the core, the fractional polarization m is in the range 
from 2.0% to 3.6%, and the average EVPA is about 110°. The 
jet component C5 shows a higher fractional polarization in the 
range form 6.0% to 6.7%, and an average EVPA -75°. At the 
last epoch 2005.19, the polarization structure appeared to dif- 
fer considerably from that observed during the previous three 
epochs. The average EVPA is about 44° and m ~ 4.8% for the 
core. The change in the polarized emission structure might be 
correlated with the total flux density variation as a result of the 
ejection of new jet components. The enhanced polarization in C5 
can be explained as the result of a compression of the magnetic 
field in this region. The compression might be caused by shocks 
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Table 8. Parameters of the hydrodynamic instability model 





Assumed parameters 




Fitted parameters 




e 


* v, y 


a 






(°) 


(°) ( mas yr l ) (mas) 


(°) 


(mas) (mas) 


(°) 


8.0 


3.0 0.19 0.0 


2.6 


0.005 0.195 


62.1 



Note: r , (po, and yo are the initial values of the distance along the per- 
pendicular jet axis, the azimuth angle, and the height along the jet axis, 
respectively; A is the spatial wavelength of the initial perturbations; 
is the velocity component along the jet axis; *F is the half-opening angle 
of the helix; 6 and a are the viewing angle and the position angle of the 
jet axis (defined as counter-clockwise from north), respectively. 

Table 9. Parameters of procession model 



Assumed parameters 

e ¥ vj 

(°) (°) (c) 


Fitted parameters 

& P tref $ jet Srot 

(°) (yr) (yr) 


8.0 3.0 0.98 


2.6 95 56 1 -1 



Note: Vj is the velocity of the jet component; the velocity vector rotates 
around the jet axis at an angle *F with a period P; s jet and s rot are two 
sign parameters used to define the jet motion direction and the helix 
rotation direction; 6 and a are the viewing angle and the position angle 
of the jet axis (defined as counter-clockwise from north), respectively. 



produced by newly ejected synchrotron plasma from the core. 
Under the compression conditions, the orientations of electron 
vector in C5 become closer to the direction of the jet flow. 

At 43 GHz, we imaged the polarization structure on a 
scale smaller than 1 mas at three epochs 2003.36, 2003.56, 
and 2004.25. The polarized emission near the core region at 
43 GHz is shown in Figure [2] At the epoch 2003.36, the polar- 
ization properties observed at 43 GHz are consistent with those 
at 15 GHz, namely EVPA^ 120° and m « 2.4%. In the second 
(2003.56) and third (2004.25) epochs, the fractional polarization 
is about 2.7% and 1.8%, respectively, but EVPA changes to 45° 
and 72°. The uncorrected variation in polarization parameters 
between 43 and 15 GHz close to the flare epoch suggests that 
the radio activity associated with the flare was embedded within 
a zone that became transparent at 43 GHz but was still optically 
thick at 15 GHz at the epochs 2003.56 and 2004.25. After ex- 
pansion, the active zone became optically thin at 15 GHz at the 
epoch 2005.19, as indicated by the significant changes in the po- 
larization properties. 



4. Helical jet structure 

Parsec- scale jets in extragalactic radio sources are often mis- 
aligned with their kilo-parsec counterparts and sh ow complex, 
in some cases "oscillating" structu re, e.g., 3C 273 (Ze nsus et al. 
]|1988[ [Conway & Murphy |1993] ). A three-dimensional helical 
pattern is often used to explain the phenomenology of oscillat- 
ing parsec- scale jets. In particular, bright knots might represent 
the locally brighter region in the jet flow caused by enhanced 
Doppler boosting effects in a helical pattern. Induced helical 
bending in jets appe ar to be common in radio galaxies and flat- 
spectrum quasars ([Conwa y & Murphy 1 1993 [[Taylor et al. |1996| 
[Pearson et al. |19981|Lobanov et al. |2Q0T^ 

Theoretical interpre tations of helical jet s include hydrody- 
namic instabilities {e.g. Hardee ||1986 1987 ) and/or precessing 



jets (e .g. |Begelman 7Blandford & Rel^[T980l |Lindfield et aT 
[T98TT ). 

In the hydrodynamic instability models, the helix might be 
triggered by small perturbations in the accretion flow or random 
perturbations, such as jet-cloud collisions at the onset of the jet. 
The initial perturbations can be amplified by Kelvin-Helmholtz 
(K-H) hydrodynamic instability. Under certain circumstances, 
the instability develops into helical modes downstream to form 
the smoothly oscillating jet seen in 2-D images. The actual evo- 
lution of the instability in the jet flow depends on the fluctuation 
properties of the initial perturbations, the dynamics of the jet 
flow, and the properties of the surrounding interstellar medium 
HHardee |1987l|2503) . 

In the precessing models, the helical jet trajectory is a su- 
perposition of ballistically moving jet knots that are for example 
ejected along periodically varying directions because of the pre- 
cession of the jet nozzle. This precession might have been caused 
by the disturbance of the je t axis by a secondary black hole in a 
binary black hole system (Begelman, Blandford & Rees |1980 ), 
or from interactions between the wobbling accretion disk and the 
spinning black hole ( [Liu &Melia |2002| ). 

In the following subsections, we analyze the jet structure of 
1156+295 in the context of both possible scenarios. 



4.1. Hydrodynamic instabilities 

A helical pattern can be produced by small perturbations in the 
jet flow by a hydrodynamic instability under certain circum- 
stances. A dispersion relation describing the behavior of the har- 
monic components of the perturbation has been derived from a 
Fourier expansion of a cylindrical jet ( [Hardee |1987| ). A particu- 
lar resonant mode can grow f rom any pert urbations in a jet flow 
via K-H instability (e. g. |Zhao et al. [T992] ) . 

To compute the oscillating jet trajectory in 1156+295 from 
sub-pc to pc scales, we introduced a kinematic helical model 
on the basis of case 2 from Steffen e t al. | ( |1995|). This model 



corresponds to the helical jet model derived by Hardee | ( [1987| ) 
based on the assumption of an isothermal and adiabatic expand- 
ing jet. For a cylindrical jet, Hardee ]( [1987| ) showed that the jet 
distortions will grow in a way primarily associated with the he- 
lical fundamental mode n =1 and therefore produce a helically 
twisted pattern on the jet's surface. The helical jet can propa- 
gate to a rather long distance, i.e., 100 times the r adius of the jet, 
prior to the disruption of the jet collimated flow ([Hardee 11979] 
1986[|1987| ). Under a relativistic approach, [Steffen et al. ||l995| ) 



calculated the a 3-D helical trajectory of the jet motion and as- 
sumed the bright knots in the helical jet are enhanced by Doppler 
boosting, where the jet plasma is likely to be either compressed 
or shocked. 

Assuming that the kinetic energy, the specific momentum 
component along the jet axis, and the opening half angle of the 
jet all remain constant, the jet trajectory of the jet m otion can be 
described by the equations of Steffe n et al. | p995 ) in a cylindri- 
cal coordinate system 



r(f) = v/tanY + r , 



0(0 = 00 + 



y(t) = v y t, 



v^tan^ 



In 



K0 



(6) 
(7) 
(8) 



where r, 0, and y are the distance along the perpendicular jet 
axis, the azimuth angle, and the height along the jet axis, respec- 
tively. We note that the notation z, the distance along the jet axis, 
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in Steffen et al. (1995) has been replaced by y, and is not to be 
confused with the cosmological redshift z. The coordinates with 
the subscript V (r , 0o> and yo) define the region where the ini- 
tial perturbations take place, v y is the velocity component along 
the jet axis, *F is the half-opening angle of the helix, and a>o is 
the initial angular velocity. 

The properties of the helix in this model are simply described 
by six parameters: four parameters (the initial angular velocity 
qjo or, the initial spatial wavelength, Aq = r , (po, and yo) 
describe the properties of the initial perturbations; and an addi- 
tional two parameters (v y and *F) describe the properties of the 
jet in which the perturbation is propagated. 

To fit the model to an observed jet trajectory, we need two 
more angular parameters to describe the projection of the 3-D 
jet trajectory on to the sky plane: the angle between the jet axis 
and the line-of- sight (or the viewing angle) and, a the position 
angle of the jet axis (defined as counter-clockwise from north). 

We first assume that the perturbation takes place in a region 
close to the central engine, so that yo can be set to zero for simpli- 
fying the calculations. The viewing angles of compact jet com- 
ponents calculated in Section 3.4 vary between 1.5° and 22.1°. 
If we inc lude published data at lower observing frequencies in 
Hong et al. (|2004|), the arithmetic mean value of the viewing an- 



We also examined the other cases considered by Steffe n et al. 
[ (T995J. In case 1, the specific angular momentum, the Lorentz 
factor, and the momentum component along the jet axis are as- 
sumed to be conserved. This case is rejected at first since no 
more than one quarter revolution can be produced for the helical 
pattern. Both cases 3 and 4 involve a magnetic field, case 4 takes 
into account the dissipation of the kinetic energy of the jet flow 
by means of collisions with the surrounding ISM and radiative 
losses, while case 3 considers the variation in the specific mo- 
mentum. However, we found that the fitted curves in the models 
of both cases 3 and 4 degenerate quickly into a straight line on 
the scale of several milliarcsecs, inconsistent with the multiple 
jet bendings observed in 1156+295 on the scales ranging from 
mas to arcsecond. 



4.2. Precessing jet 

Alternatively, the observed helical jet morphology might repre- 
sent a superposition of ballistic jet co mponents ejected seq uen- 
tially from a precessing jet nozzle (e.g. |Stirling et al. |2003) . On 
the basis of the assumption that the jet components are ejected 
from a precessing jet nozzle and move along ballistic paths, we 
fit the bright knots observed in the jet with a simple precessing 



gles of compact jet components is 8.0° (+4.9). We take the mean 
value as an appropriate estimate of the viewing angle of the helix 
axis 6. The half-opening angle *F of 1 156+295 is estimated to be 
~ 3° (determined from the maxima range of the position angles 
and viewing angles of bright jet components). 



Furthermore, we derived the jet velocity Vj in Section [3T4 



Since both the opening angle of the helix cone and the view 
ing angle of the jet are small, the velocity along the jet axis 
is approximately equal to the jet velocity, « vj ~ 0.98 c. 
Taking into account of the effect in special relativity, the ap- 
parent jet velocity is boosted by a factor of 



1 - v y cos 6/c 



mode l (Hjellming & Johston |198T||Baryshev |1983t[AneFaT 
|20T0| ). 

In this model, the jet knot moves outward from the central 
object with a constant velocity vj. The velocity vector Vj rotates 
around the jet axis at an angle *F and a period P. Each jet com- 
ponent travels along a ballistic trajectory, and the projection of 
the superposition of their ballistic trajectories constitutes the ob- 
served helical jet trajectory on the sky plane. The angle between 
the jet axis and the line of sight is and the observed jet axis lies 
at a position angle a. 

The velocity o f a given jet can be described as in Hjellming 
|& Johston I ( [1951) with 



4.7, i.e., Vy P = v y 



sin# 



I -v y cos 6/c 



0.19 mas yr" 1 . 

The remaining four parameters (a, r , (po, A) are the free pa- 
rameters to be determined by fitting the model to the observed jet 
trajectory. As illustrated in Figure [9] the observed jet trajectory 
can be well fitted with the isothermal helical jet model on wide 
scales ranging from 0.08 mas (0.6 pc) to 35 mas (250 pc). The 
symbols represent the jet components detected from the 15 and 
43 GHz VLBI data from 2003 to 2005 presented in this paper, 
as well as the ones detected at 1.6, 22, and 86 GHz from other 



4.6 c, or a proper motion of v x = s jet v[sm *F cos £l(t - t ref ) cos 6 + cos *P sin 0], (9) 

v y = Sj et v sin *F sin Q(t - t re f), (10) 



*jet 



v[cos *F cos 6 - sin *F cos £l(t - t re f) sin 6], 



(11) 



papers ( |Hong et al. |2004| |Jorstad et al. |2001[ |Lee et al. |2008] ) . 

The best-fit model parameters are given in Table [8] The po- 
sition angle a « 2.6° in the source frame is inferred from our 
model fitting, and is consistent with the observed jet structure 
that extends northward from the core. 

Our model fitting suggests that there is a region with r ~ 
0.005 mas (or -0.04 pc) in which the jet flow might be perturbed 
substantially. This compact region in the vicinity of the central 
engine might be close to or associated with the jet nozzle. This 
radius (0.04 pc) is about 1000 times of the Schwarzschild radius 
(R s « 4 x 10" 5 pc) of a supermassive black hole (SMB H) of 
4.3xl0 8 M o which has been inferred to exist in 1156+295 (|Pian 



where Sj et is used to define the jet motion direction and Sj et = +1 
corresponds to moving towards the observer, t re f is the reference 
time, D is the rotating angular velocity of the jet velocity vector, 
Q = s rot 2n/P, s rot is used to define helix rotation direction, and 
s rot = +1 corresponds to counterclockwise. 

The transverse velocity of the jet can be decomposed 
into two orthogonal components of right ascension (a) and 
declination^) 



vs = cos avy - sin av z . 



(12) 
(13) 



Taking the relativistic aberration into account, the proper 
motion of a given component ejected at time t e j ect on the sky 
plane can be written as: 



et al. ||2005| ). It appears to be much greater than the inner edge 
of the accretion disk radius of 3 R s , equivalent to the radius of 
the inner most stable orbit, but comparable to the typical size of 
the BLR. Perturbations in the jet flow might be imposed by the 
interactions between the jet plasma and BLR clouds. 



v a (t -t eject )/ COS S 
d(l-v x /c) ' 

Vd(t - t e ject) 

d(l-v x /c)' 



(14) 
(15) 
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where d is the distance from the observer to the source. 

A jet velocity Vj > 0.98c is derived in Section 3.4 in the 
source-rest frame. The jet half-opening angle *F of 3° and view- 
ing angle of 6 = 8° are the sa me a s the values adopted in our 
K-H instability model (Section pTT). Since the jet is moving to- 



ward us, Sj e t is set as +1. Adopting an eyeball approach, we fit 
the precessing model to the jet structure observed on the parsec 
scale by adjusting the remaining free parameters (position angle 
of the helix a, precessing period P, reference time t re f, and the 
sign of rotation s rot ). 

The fitted parameters are given in Table [9] As illustrated in 
Figure [T0| on the scale of tens of mas, the precessing jet model 
provides a reasonable fit to the observed data but the inner 2- 
mas jet structure appears to be poorly fitted. The structure on 
the scale > 35 mas is also difficult to fit with the precessing jet 
model. The kinematic model described here represents an ideal 
case, in which the jet is continuously ejected at a constant veloc- 
ity and the jet ejection direction varies periodically. However, the 
actual jet precession in a powerful radio source might be more 
complex. For example, the period of the jet ejecting direction 
may be altered by intermittent radio outbursts and jet ejections 
( Reynolds & Begelman 1997| ). The jet may interact with the 
BLR and/or NLR clouds to generate stationary shocks. Abrupt 
changes in the jet direction and violent jet-ISM interactions have 
been o bserved in som e compact-steep-spectrum sources such as 
3C 48 flAn et al. 12010] ) and 4C 41.17 ( [Gurvits et al.TT997) . In 
addition, the intrinsically small jet bending in the inner 2 mas 
might be enhanced by projection effects. There is also a possi- 
bility that hydrodynamical instabilities of the plasma at the onset 
of the jet may modulate the fine structure of the whole curved jet. 

A precessing jet from AGNs is a characteristic phenomenon 
in a binary black hole (BBH) system (Begelman, Blandfo rd &| 
Rees [1980 ). However, the orbital period of the BBHs remains 
uncertain because the relationship between the precessing pe- 
riod and orbital period is unclear. Using the ratio of the preces- 
sion period to the orb ital period in X-ray binaries as an analogy 
(Marg on et al. 11980] ), the precessing period fitted for 1156+295 
gives a upper limit to the orbital period of the BBHs of about 95 
years. 

Assuming Keplerian motion, the derived semi-major axis, or 
the largest separation between the two black holes, is <0.08 pc 
for a primary black hole with a mass of 4.3xlO 8 M . This size 
scale is one tenth of the current VLB I resolution. Moreover, the 
real orbital period is probably a few percent of the precessing 
period, thus the orbital separation might b e a few xO.Ol pc. 

Some theoretical calculations (e.g. |Yu |2002] ) show that the 
orbital parameters of a survival BBH depend on the galactic ve- 
locity dispersions and blackhole mass ratios. An orbital period 
of a few tens of years suggests a mass ratio of the secondary 
to primary black holes of between 0.01 and 0.05. Combining the 
close separation of the binaries and the relatively low mass of the 
secondary BH, the radio emission from the secondary BH might 
be very weak and below our current detectable limit in VLBI. 
Evidence of possible BBHs in 1156+295 may be represented 
by possible periodic variabilities in the radio, optical, and X- 
ray light curves, and the expected apparent variability timescale 
would be suppressed by relativistic aberration resulting in an or- 
bital period much shorter than the intrinsic value ( |Yu |2 002). 



4.3. Comparison of the two models 

In the hydrodynamic instability model, the jet trajectory on 
scales from sub-mas to several tens of mas and even larger can be 
naturally explained as the pattern produced by the fundamental 



mode (n = 1) or a helical mode in K-H instability model with an 
initial characteristic wavelength of Aq = 0.2 mas (1.5 pc). On the 
basis of our model fitting to the data, the helical mode appears to 
be excited at the jet base on a scale of 0.005 mas or 10 3 R S , within 
the scope of the BLR. The emission from the jet knots moving 
towards the direction of observers is Doppler boosted. In com- 
parison, the precessing jet model alone appears to have difficulty 
in fitting the jet trajectory on the observed scales ranging from 
sub-mas, several tens of mas, up to arcsec. In particular, the mul- 
tiple bends observed on the mas scale cannot be explained solely 
with the precessing jet model. With additional mechanisms such 
as interactions of jet plasma with the interstellar medium, we 
may be able to explain the considerable deviation of the observed 
jet trajectory from that predicted in the precession jet model. 
Furthermore, the jet trajectory observed beyond 35 mas imposes 
a further limit on the precessing jet model using a single rotat- 
ing angular velocity (Q) of the precessing jet nozzle. To fit the 
precessing jet model to the data observed beyond 35 mas sug- 

2tt 

gests that an additional rotating period (P = — ), which is much 

longer than 95 years derived from the inner jet, is required. For 
1156+295, despite the deceleration, the jet plasma on pc scales 
beyond 35 mas has been ejected from the jet nozzle for at most 
a few hundred years. It would be incredible if the period of a 
process associated with the nuclear dynamics had changed con- 
siderably and abruptly in such a short period. Finally, based on 
our VLA observations ( Hon g~eTal. |2004| ), the jet trajectory be- 
comes a straight line on the kpc scale (or >0.1 arcsec). In con- 
trast, the precessing jet model predicts that the curvature of the 
jet trajectory projected on the sky becomes large unless ballis- 
tic jet components ejected from the precessing nozzle in a wide 
range of directions had been lined up, on the way travelling to 
the outer region on the kpc scale, by a tunnel in the ISM or by a 
poloidal magnetic field in 1156+295. 

In short, although the results from our model fitting support 
the idea that the K-H instability is responsible for the observed 
helical jet pattern, we cannot fully rule out the possibility that 
the oscillating jet structure observed between 1 mas (10 pc) and 
40 mas (300 pc) is driven by a precessing jet nozzle with a ro- 
tating period of 95 yr. Additional astrophysical processes (e.g. 
interaction between the jet plasma with the the BLR and/or NLR 
clouds) should be taken into account to explain the complex 
bending jet within 10 pc, and recollimation of the jet beyond 
300 pc, if the jet precessing model is to be successfully applied 
to 1156+295. 



5. Summary 

We observed 1156+295 with the VLBA at the wavelengths of 
86, 43, and 15 GHz with angular resolutions up to 0.08 mas 
(0.5 pc) at four discrete epochs during 2003 to 2005. Based on 
the analysis of the new VLBA data from the multiple-epoch and 
multiple-frequency observations, we have presented new results 
on the morphology, physical properties, and kinematics of the 
jet observed on the scales from sub-pc to pc from the blazar. We 
have also discussed the astrophysical implications for the helical 
pattern of the jet emission between sub-pc and kpc by fitting both 
the K-H instability model and the jet-nozzle precessing model to 
the observed data. 

The core-jet structure was detected in 1156+295 at 43 and 
15 GHz at angular resolutions from 0.2 mas to 0.6 mas. Six 
jet components were identified. Among the six, the innermost 
three were detected for the first time. The apparent transverse 
velocities of the six jet components range between 3.6 c and 
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1 1.6 c. At 86 GHz, during the first epoch the source was resolved 
into a core-jet structure on the image, during the last epoch the 
source seems to be resolved into a double-point source in the uv- 
domain, and during other epochs only a point-like source was 
detected on both the image and the wv-domain. Ultra-high appar- 
ent brightness temperatures of > 10 12 /<r were inferred suggest- 
ing that highly relativistic jet plasma ejected from the nucleus 
moves in the directions close to the line-of-sight of observers. 
After correcting for both the Doppler beaming and cosmolog- 
ical redshift effects, the intrinsic brightness temperature of the 
nuclear radio emission is below the inverse Compton catastro- 
phe limit. During the observing period from 2003 to 2005, a ra- 
dio flare was present in the data, suggesting that a radio activity 
of the blazar originated inside. Linear polarization was detected 
at 43 and 15 GHz in both the core and inner jet. During the qui- 
escent epochs, the EVPA of the jet polarized emission differed 
by ~ 35° from that of the core. No significant change in EVPA 
of the polarized emission from the core between 2 cm and 7 mm 
was observed during the quiescent epochs. At the epochs after 
the radio flare, EVPA observed at 7 mm from the core shows a 
significant change from that of the quiescent epoch. 

By model fitting a helical fundamental mode of K-H insta- 
bility and a precessing jet-nozzle, we have found that the K- 
H instability model successfully fits the observed jet trajectory 
on the scale ranging from sub-pc to kpc. Our results suggest 
that the fundamental helical mode with an initial characteristic 
wavelength of 0.2 pc appears to be excited at the jet base on 
a scale of 0.005 pc or 10 3 #„ a typical BLR size for a SMBH 
of 4.3xl0 8 M©. The precessing jet model only fits the jet tra- 
jectory observed on a scale ranging from 10 pc to 300 pc and 
requires additional astrophysical processes to explain the com- 
plex jet bending observed in the inner region (sub-pc to 10 pc) 
and the re-collimation of the large scale jet (> 300 pc). 
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Fig. 1. The four epochs of VLBA images of 1156+295 at 15 GHz. The contour maps represent Stokes I images. The pseudo-color 
maps overlaid on contours represent the distribution of fractional polarization in which linear polarization intensity above 4<x is 
shown. The length of the short bars represents the strength of the polarized emission, 1 mas~1.0xl0~ 2 Jy beam -1 . The orientation 
of the short bars represent the apparent electric vector position angle (EVPA) without a rotation measure correction. The detailed 
parameters of the images are listed in Tables [2] and [7] 
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Fig. 2. The four epochs of VLBA images of 1 156+295 at 43 GHz. The contour and pseudo-color maps are similar to those in Figure 
[I] The length of the short bars represents the strength of polarized emission, 1 mas~4.0xl0~ 2 Jy beam -1 . The orientation of the 
short bars present the apparent EVPA without rotation measure correction. The detailed parameters of the images are listed in Tables 
ElandEl 
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Fig. 3. The fitted models are overlaid on the clean images of 15 and 43 GHz on epoch 2003.36 to indicate their positions and intrinsic 
size. 
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Fig. 4. Estimating the flux density (S f) and size (6d) of the source using the brightness distribution models introduced by Pearson 
(1995 ): for the first three epochs, an extended source model could be used, whereas for the last epoch, a double-point source model 
with an intensity ratio of 3:1 and an angular separation of 0.16 mas were used 
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Fig. 5. The four epochs VLBA images of 1 156+295 at 86 GHz. 
The detailed parameters of the images are listed in Table [6] 
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Fig. 7. The proper motions of the outer jet components C2, 
C3, and C4 are 0.49+0.03 mas yr~\ 0.41+0.08 mas yr" 1 , and 
0.18+0.03 mas yr" 1 , which correspond to 11.6+0.7 c, 9.7+1.9 c, 
and 4.3+0.7 c, respectively. 
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Fig. 8. The proper motions of inner jet components C5, C6, 
and C7 are 0.15+ 0.02 mas yr -1 , 0.19+0.02 mas yr -1 , and 
0.16+0.04 mas yr" 1 , which correspond to 3.6+0.5 c, 4.5+0.5 
and 3.8+0.9 c, respectively. 
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Fig. 6. A spectral plot of the core component at the four epochs. 
The error bars show errors in the measurement of the flux density 
of the core component. 
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Fig. 9. The best fitting of a model with the fundamental helical mode in K-H instability (red curve) to the data (symbols). The color 
of symbols represents the data taken at different frequencies: blue for 15 GHz or lower frequencies, magenta for 22 GHz, red for 
43 GHz, and green for 86 GHz. In addition, different types of symbol mark different jet components: star for CO, cross for CI, open 
circle for C2, open diamond for C3, open triangle for C4, open pentagram for C5, open hexagram for C6, and open square for C7. 
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W. Zhao: Sub-pc structure of blazar 1156+295 
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Fig. 10. The best fitting of the precessing-jet model (red curve) to the same observed data as described in Figure [9] The color of 
symbols represents the data taken at different frequencies: blue for 15 GHz or lower frequencies, magenta for 22 GHz, red for 
43 GHz, and green for 86 GHz. In addition, different types of symbol mark different jet components: star for CO, cross for CI, open 
circle for C2, open diamond for C3, open triangle for C4, open pentagram for C5, open hexagram for C6, and open square for CI. 
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